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1 INTRODUCTION 



At low energy /i ~ m T an expansion parameter of perturbation theory the strong coupling 
constant a s = a s (m T ) ~ 0.35 is rather large (e.g. Q) and, therefore, for a generic QCD 
observable a (without Born term) 

a ~ a s (l + aia s + a 2 a 2 s + . . .) (1) 

the series of perturbation theory approximation converges badly the last term being about 
10% of the leading one. 

To improve predictions that is required by experimental data for some observables at 
present, higher order corrections have to be included. This, however, is technically difficult 
due to the necessity of computing many-loop integrals representing Feynman diagrams in 
high orders of perturbation theory. There is a little hope to obtain next terms in (|l]) 
(beyond three-four loops) for many processes in realistic models. 

It should be stressed that phenomenologically (leaving apart general considerations 
of possible power corrections that are necessary for two point correlators and lead to 
resonances) there is no much intrinsic reason to go beyond perturbation theory: wild 
asymptotic behavior is not seen yet. Also the freedom of the choice of the renormalization 
schemes allows one to render the series convergent for a given observable (or some set of 
observables) at least in an heuristic sense that further known terms decrease [||. There is 
no strict indication that perturbation theory is broken though the accuracy it can provide 
in a number of cases is not sufficient for confronting predictions with experimental data. 
Thus, at the level of phenomenology just needs of precision require an improvement of 
perturbation theory predictions and because further terms are not available going beyond 
perturbation theory in different ways is now widely discussed Q. 

Before adding genuine nonperturbative terms which is not obvious in cases when Wil- 
son operator product expansion is not directly applicable one tries to go beyond the finite 
order perturbation theory by a resummation of a particular subset of terms that can be 
explicitly generated. The simplest one is due to running of the coupling constant. 

This resummation is ambiguous to a great extent in particular it can change the 
analytic properties that exist in any finite order of perturbation theory and are established 
on a general ground of quantum field theory and even can make them wrong i.e. 
resummed quantities can not satisfy some general requirements that leads to necessity of 
interpretation of the results. The ambiguities that are produced by the resummation and 
the change of analytic properties are analyzed in some details. 
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2 AMBIGUITY FOR THE r LEPTON WIDTH IN MS 
SCHEME 



The r lepton width became a real laboratory for investigation of properties and numerical 
validity of low energy perturbation theory || . 

The spectral density R(s) for a two point correlator 

U(x) = <0|Tj(x)j(0)|0) 

where j(x) is a weak charged current of light quarks generates Adler's function 



D(Q 2 ) = Q 2 I j^2>Q 2 = -* 2 (2) 



that can be calculated in Euclidean domain in terms of perturbation theory series in the 
coupling constant a s (n) 

D(Q 2 ) = a( M ) + aififiPo In ^ + c) + . . . 

From the last expression the spectral density can be found in any finite order of pertur- 
bation theory in the form 

„2 

R(s) =a(ji) +a{n) 2 (p In ^ + c) + . . . (3) 

s 

The r lepton decay width are given by 

r T = / R(s)W(s/m 2 T )ds/m 2 , (4) 
Jo 

where 

W(x) = (l-x) 2 (l + 2x). 
In a finite order of perturbation theory the expression for the r lepton width has the form 

r> = a(u) + a(//) 2 (/3 ln — ^ + c) + . . . 

Above formulas are given for normalization only and /?o is the first coefficient of the (5 
function in QCD. Within the formal perturbation theory in the strong coupling constant 
a s every term (ln n -^) has correct analytic properties in g 2 : a cut along the positive 
semiaxes. Using formulas @ and (Q) with known renormalization group properties of R(s) 
and r r (/i-independence) one can go beyond the predictions of finite order perturbation 
theory and perform a partial resummation in many different ways. First it can be done 
directly on the cut ||. Because the spectral density itself is nonintegrable at low energy 
after using the renormalization group improvement on the cut (fi 2 = s) 

i 

R(s) ~ a s (s) 



PoHs/A 2 ) 
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one can take the discontinuity after RG summation in Euclidean domain to obtain (in the 
leading order) 

1 7T 2 

R(s) = — arctan ira(s) = a(s) a(s) 3 + . . . (5) 

7T 3 

Now the integral in (||) can be done that provides an improvement of PT prediction 
through resummation on the cut. This technique works up to third order of PT that is 
available now, i.e. commutation of RG summation and taking the discontinuity makes the 
integral in (||) regular. 

At third order of perturbation theory, however, there appears another way of definition 
of the quantity in question (||) that is connected with the change of renormalization scheme 
that reduces to a redefinition of the coupling constant. In the effective charge scheme ]?]] 
an effective j3 T {a T ) has a zero that leads to an IR fixed point || and integrals (||) can be 
explicitly done in this scheme ||. This possibility however depends crucially on the order 
of PT and is absent in the second order. It is unknown whether it persists in the fourth 
order. So for this technique the natural requirement that the method of resummation is 
stable in every order of PT is not fulfilled. 

There are also other possibilities of resummation, for instance, different kinds of opti- 
mization || . 

Another recipe (more perturbative because it is formulated in the complex plane and 
not on the physical cut) is to use the analytic properties given by (|2|) and define |11] 

/ R(s)ds = — f U(z)dz. 
J 2m Jc 

Cauchy theorem requires no singularities inside the contour so if 

D{z) ~ D(z) PT ~ a(z) = / 

/5 ln(-z/A 2 ) 

then D(z) PT has wrong analytic properties and there is a difference (nonperturbative) 
which is proportional to A = m T exp(— 2-7r/9a s (m T )) with the result of direct summation 
on the cut. 

One should stress that the modified minimal subtraction scheme is always used for the 
definition of the charge and the change of an observable within PT is formally of higher 
order in the coupling constant. The real problem is that this difference is large enough to 
be caught by experiment. Then the theoretical predictions can differ by the amount that 
depends on the procedure used in computation and is not negligible. It is unclear how to 
single out the best numerical value. There is here even more ambiguity than the simple 
freedom in the choice of the renormalization scheme. 

Thus, aMs( m r) depends rather strongly on the resummation procedure that should 
be explicitly explained when precise comparison of different predictions is made. 

Having this ambiguity in mind and noting that (x^(m T ) by itself can not be measured 
because it is unphysical quantity we next consider more strict test of pQCD that involves 
only observables and therefore is free of the renormalization scheme ambiguity. Still an 
ambiguity due to resummation contrary to the finite order analysis is present. 
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3 A TEST OF pQCD FOR DIRECTLY MEASURED 
QUANTITIES 

For the analysis the moments of e + e~ annihilation and r T are chosen because [12, 13] 

• these observables are generated by the same Green's function in pQCD (m q = 0) 
that reduces unknown possible nonperturbative effects 

• the integration scale for the moments can be adjusted in such a way to avoid the 
renormalization group evolution that would require the use of the (5 function and 
introduce further uncertainties 

• both moments and r T can be directly measured with high precision that allows to 
pin down the theoretical difference that is parametrically of the next (fifth) order in 
the coupling constant and is fairy small. 



Notations for further analysis are as follows [13]. The whole spectral density 



R{s) = 2(1 + -r(s)) 



is defined through the reduced one 

9 

r(s) = —a + ... 
Air 

that determines the reduced Adler's function 

r(s)ds 



Jo {s + Q 2 ) 2 
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+ M- +k n~ +M- +••• 

The moments of e + e" annihilation rate are defined by 

f m r ds ( s \ n . . 
r n = (n + 1) / — 2 — r{s) 
Jo m 2 \m4J 

r T = 2r - 2r 2 + r 3 . 

Coefficients ki summarize all information from perturbation theory and are the only 
ingredient for testing the theory. We factor out all known nonperturbative corrections due 
to condensates. The technique of resummation based on integration along the contour in 
the complex plane is adopted |ll]] • Because r T and r n are renormalization group invariant it 
is convenient to use renormalization group invariant approach from the the very beginning. 
Introduce 



d T = d(m^) = m 2 







r(s)ds 
(s + m" 2 .) 2 
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for which the RG equation is 

z^-d = -d 2 (l + pid + p 2 d 2 + p 3 d 3 + ...) 
dz 

where pi are renormalization scheme invariants, p\ = 0.79, p 2 = 1.035, P3 = 2k% — 2.97953 
(MS parameterization) with recently computed coefficient of the j3 function J[4|]. 

Moments are defined as integrals along the contour in the complex plane and therefore 
unphysical singularity is included. If 

p(z) : - Z ^P( Z ) = d (z) 



then 

n+1 



x n p(m T x)dx. 

x\=l 



2ni 

The machinery of computing consists now in finding functions f n (-) and g(.) and then 
inverting the function g(.) to represent moments through the only input parameter r T 

r n = fn(d T ), r T = g(d T ), d T =g~ 1 (r T ), 

r n = /nfo -1 (»Y)] = {fn® <? _1 )K)- 
Analytic properties of the function f n ® 9~~ 1 (-) in r r determine a radius of convergence for 



our series [|12[. We could not determine these properties completely. However the simpler 
question - analytic properties of f n (x) with respect to x - can be answered completely (in 
MS scheme without k 3 ). To the leading order the explicit formula is 



that gives x<l/7rora<|^ [|D| . 

Generalization to higher orders is straightforward, one has to find a singularity of the 
solution of the renormalization group equation ||] 

z-^d(z) = p(d(z)), d(m 2 T ) = d T . 

The result is < §0.744, < |o.697, < f 0.674 @. 

The actual value of a (3) = 0.3540 obtained from r e T xp = 0.487 ± 0.011 @ lies outside 
convergence regions. The behavior of the decay rate in higher orders of perturbation 
theory is presented in Table [l] If the resummation does catch a dominant behavior 
then next several orders of perturbation theory do not improve much and do not show 
wild asymptotic behavior either. It forms a kind of dead zone where no qualitatively new 
behavior starts. 

Results for the moments are presented in Table |2[ The label "PT" stands for finite 



order perturbation theory results [12], while "^2" and ll p^ n columns contain resummed 



results in third and fourth order approximation for the (3 function [13]. "Exp" gives 
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n 


r T 


n 


r T 


n 


7> 


1 


0.2535 


5 


0.4993 


9 


0.7400 


2 


0.4021 


6 


0.4622 


10 


1.5577 


3 


0.4870 


7 


0.4402 


11 


4.059 


4 


0.5153 


8 


0.4894 


12 


11.905 



Table 1: Higher order perturbative predictions for the reduced part of the semileptonic r 

(3) 

decay width using a T = 0.2535 





PT 


P2 


P3 


exp 


Ro 


2.28 ±0.05 


2.334 


2.338 


2.15 


Ri 


2.14 ±0.08 


2.219 


2.216 


2.06 


R2 


2.12 ±0.12 


2.228 


2.230 


2.00 


er- 


trunca- 


~ 0.07 


~ 0.07 


- 10% 


ror 


tion 


input 


input 


sys 



Table 2: Predictions for moments of e + e annihilation rate through r lepton width 



experimental number after direct integration of experimental data taken from the compi- 
lation [^] . Last line indicates main source of the uncertainties which are due to truncation 
of the perturbation theory series, an input value of r^ xp and systematic errors of the data. 

Note that though we take into account the running of the coupling constant that is 
sensitive to an IR domain, due to a particular way of analytic continuation there is no 
direct problem of infrared renormalons J18[ |. 

At a given order of the (3 function theoretical predictions with resummation contain 
no error but due to input parameters (r T in this case) . Errors of summation method itself 
cannot be strictly given however. As an estimate we use the change of predictions when 
passing from one order to the next one, i.e. p2 — > P3- To find p% one needs the numerical 
value of £?3. For the popular estimate of ks [[0], IS] one has p% = 1.36. After the analysis 



done in [13] we propose a slightly different number p^ = 2.0 ± 0.5 that we have used for 
our predictions (recall that in our normalization p\ = 0.79, p% = 1.035). 

The only quantity that is given at the finite order of perturbation theory in our ap- 
proach is the (3 function and one has to check its validity at least heuristically in terms of 
decreasing with order. The worst pattern of convergence for the (5 function in the course 
of the analysis is given by (numerically) 

(3{a) ~ 1 + 0.284 + 0.134 + 0.0467p 3 + . . . 
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at the point a = 0.36. In all other points on the contour the convergence is better. One 
can see that the numerical value of p% is rather important for application of discussed 
technique even though the change of results is small. 



4 CONCLUSIONS 



As conclusions to my talk I summarize results of investigation done in 0, O, 13]: 



the coupling constant extracted from different processes even in the same renor- 
malization scheme (for instance, MS) requires for precise comparison an explicit 
mentioning of the resummation procedure 

the recipe based on analytic continuation is rather stable against inclusion of higher 
order corrections to the [3 function 

moments of e + e~ annihilation at /i ~ m T are estimated to be smaller than required 
by pQCD with resummation of effects of running of the coupling constant. 
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